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We study the 7-ray emission properties of pulsars by using a new self- 
consistent outer gap model. The outer gap can exist in pulsars with age over 
million years old if the effect of magnetic inclination angle as well as the average 



i properties of the outer gap are considered. The mature 7-ray pulsars, whose 

^ I ages are between 0.3 to 3 million years old, are able to move up to high galactic 

2 I latitude. Moreover, their 7-ray luminosity are weaker and their spectra are softer 

Tjj- I than those of younger 7-ray pulsars in the galactic plane significantly. We use 



a Monte Carlo method to simulate the statistical properties of 7-ray pulsars in 
the Galaxy as well as in the Gould Belt. We find that 7-ray pulsars located at 
I 6 I < 5° and located at | 6 | > 5° have very different properties. High galactic lat- 



O 

H i itude 7-ray pulsars are dominated by mature pulsars with longer periods, weaker 



O ■ 

^ I fluxes and softer spectra. If the pulsar birth rate in the Galaxy and the Gould 

Belt are ~ lO^^yr^^ and ~ 2 x lO^^yr^^ respectively, there are 42 and 35 radio- 
quiet 7-ray pulsars for I 6 |< 5° and | 6 |> 5° respectively. Radio-quiet 7-ray 



' pulsars from the Gould Belt are 2 and 13 for | 6 |< 5° and | 6 |> 5° respectively. 

We suggest that a good fraction of unidentified EGRET 7-ray sources may be 
these radio-quiet 7-ray pulsars. Furthermore 7-ray pulsars located at | 6 |> 5° 
satisfies oc L^^ whereas oc Lf^ for 7-ray pulsars in the galactic plane, where 
P ~ 0.6 and 6 ~ 0.3 respectively. 

Subject headings: gamma-rays: theory - pulsars: general - stars: neutron - stars: 
statistics 



1. Introduction 



There are 170 unidentified 7-ray sources in the third EGRET catalog (Hartman et al. 
1999), where ~50 sources close to the Galactic plane with \b\ < 5° and ~70 sources in 
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the medium latitudes with |6| between 5° and 30°. For those unidentified 7-ray sources in 
the Galactic plane, many of them are associated with Wolf-Rayet and Of stars, SNRs and 
OB stars(Montmcrlc 1979; Kaaret & Cottam 1996; Yadigaroglu & Romani 1997, Romero 
et al. 1999). Most of these objects are considered as pulsar tracers, therefore it is natural 
to suggest that these low latitude sources may be Geminga-like pulsars, which are radio- 
quiet pulsars (Yadigaroglu & Romani 1995; Cheng & Zhang 1998; Zhang, Zhang & Cheng 
2000). Since 7-ray pulsars are known to be steady 7-ray emitters, it has been suggested 
that the variability should be a good indicator to identify the real pulsar candidates from 
the unidentified EGRET sources (Romero, Combi & Colomb 1994; Mclaughhn et al. 1996; 
Zhang et al. 2000; Torres et al. 2001; Nolan et al. 2003). Gehrels et al. (2000) define a 
class of sources as steady: a source is steady if the most significant detection of a source in 
3EG catalogue is for a timescale of years and if that particular flux is within 3a of the flux 
calculated for the full data set. From this classiflcation, 48 unidentifled sources at |6| < 5° 
and 72 unidentified sources at |6| > 5° arc steady sources respectively. On the other hand, 
Grenier (2001) suggests another definition called persistent sources, which are those detected 
sources with a significance \/TS > 4 for every observation at |6| > 2.5°. If we assume that 
the sources with a significance \/TS > 5 at |&| < 2.5° are persistent sources, there are 40 
(45) persistent unidentified EGRET sources at |6| < 5° (|6| > 5°). 

However, the spectral properties of medium latitude sources are significantly softer, 
fainter and have a steeper logN-logS function than those at low latitudes (Gehrels et al. 
2000). It has been suggested that they are associated with recent supernovae in the nearby 
Gould Belt (Grenier 1997; Gehrels et al. 2000; Grenier 2000). Their natures remain as 
mystery. Harding and Zhang (2001) used the polar cap models (Daugherty & Harding 1996; 
Harding & Muslimov 1998) to investigate whether 7-ray pulsars viewed at a large angle 
to the neutron star magnetic pole could contribute to unidentified EGRET sources in the 
medium latitudes associated with the Gould Belt. They suggest that the off-beam 7-rays 
come from high-altitude curvature emission of primary particles can radiate over a large solid 
angle and have a much softer spectrum than that of the main beams, and at least some of 
radio-quiet Gould Belt sources detected by EGRET could be such off-beam 7-ray pulsars. 

Recently, the brightest of unidentified EGRET sources in the medium latitude, 3EG 
J1835-I-5918, is strongly suggested as the second Geminga-like pulsar ( radio-quiet pulsar) 
by multi-wavelength observations including Chandra, HST and Jodrell Bank (Mirabal et al. 
2000; Mirabal & Halpern 2001; Halpern et al. 2002). The X-ray spectrum can be described 
by two components: a soft thermal X-ray spectrum with a characteristic temperature (Tinf 
3 X lO^K and a power law hard tail with a photon index 7 Ri2, which closely resembles to the 
X-ray spectrum of Geminga. The repeated radio observations at Jodrell Bank did not show 
any periodicity. The X-ray data suggests that the distance is between 250-800pc, which is 
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consistent with the distance to the Gould Belt. 

In this paper, we use the revised outer gap models (Zhang et al. 2004) to investigate 
emission properties of 7-ray pulsars in the Galaxy as well as in the Gould Belt. This revised 
outer gap model is based on the original outer gap models (Cheng, Ho & Ruderman 1986a, 
hereafter CHR I, 1986b; Zhang Sz Cheng 1997) but takes into account the effect of the 
inclination angle (a), which is the angle between the magnetic axis and the rotation axis, 
in determining the size of the outer gap, which is defined as the ratio between the dimension 
of the gap perpendicular to the magnetic field and the light cylinder radius. The revised 
model also takes into account the fact that the typical radiation region of the outer gap 
is not necessary at half of the light cylinder instead it should be better represented by an 
appropriate average over the entire outer gap. This effect is particular important for old 
pulsars. If the outer gap is only represented at half of the light cylinder, then the outer gap 
is assumed to be turn off when the gap size at this region is larger than unity. This new 
model has taken the entire active region of the outer gap into account. As long as the gap 
size is less than unity in some parts of magnctosphere, the outer gap still exist. This effect 
allows some pulsars with appropriate combination of a, P and B, can maintain their outer 
gaps until a few million years old. These pulsars are able to move up to high galactic latitude 
and their ages make them weak 7-ray sources. Furthermore, we will show that these pulsars 
will emit softer spectra and their relation between 7-ray luminosity and spin-down power 
differs from that of younger pulsars located in the galactic plane. We organize the paper as 
follows. In section 2, we review the new outer gap models. In section 3, we study the 7-ray 
emission properties of pulsars by using the new self-consistent outer gap model. In section 
4, we describe a Monte Carlo simulation method to determine the statistical distributions of 
radio-quiet and radio-loud 7-ray pulsars. In section 5, we summarize the simulation results 
and discuss their implications. Finally, a brief conclusion is given in section 6. 



2. Outer Gap Models 

Zhang & Cheng (1997) have proposed a self-consistent mechanism to describe the high- 
energy radiation from the rotation-powered pulsars. In this model, the radiation mechanism 
of relativistic charged particles from a thick outcrmagnetospheric accelerator (outer gap) is 
synchro-curvature radiation (Cheng & Zhang 1996) and the characteristic energy of high- 
energy photon emitted from the outer gap is determined by the pulsar global parameters, i.e. 
rotation period P and the dipolar magnetic field B, as well as the fractional size of the outer 
gap (/), which is a ratio between the mean vertical separation of the outer gap boundaries 
in the plane of the rotation axis and the magnetic axis to the light cylinder radius, and is 
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given by 

/ N -13/8 

E,{f) ^ 5.0 X 10' f/'B%'p-y' (^£J eV (1) 

where B12 is the dipolar magnetic field in unites of lO^^G, Rl is the light cylinder radius 
and r is the distance to the neutron star. The 7-ray spectrum becomes exponentially decay 
beyond E^{f). The fractional size of the outer gap determines the total 7-ray luminosity 
from the outer gap (Zhang & Cheng 1997) and is given by 

-^7 — f^Lsd (2) 

where Lg^ = 3.8 x 10^^Bf2P~*ergs^^ is the spin-down power of pulsar. CHRI has shown 
that the inner boundary of the static outer gap begins at the null charge surface(r2 -6 = 0). 
Recently Hirotani & Shibata (2001) has shown that if there is injected current from the 
inner boundary, the position of the inner boundary can shift either toward the star or even 
close to the light cylinder depending on the sign and the magnitude of the injected current. 
However, it is not clear what causes this injected current. In this paper we shall assume that 
the outer gap begins at the null charge surface. The radial distance from the null surface 
surface to the star is r.^, which is a function of the magnetic inclination angle a and is given 
^ ~ sin™ sin^V"-") ' ^'^^^^ polox angle between of rj„ determined by tan^j„ = 

^(3 tan a + V9 tan^ a + 8) and 9^ is the polar angle to the position where the first open field 
line intercepts the light cylinder determined by tan^^c = — ^^^^^ (1 + (1 + 8 tan^ Q/g)-*-/^). The 
fractional size of the outer gap is limited by the pair production between the soft thermal 
X-rays with characteristic energy Ex from the stellar surface and the high-energy photons 
with energy E^{fo) emitted from the outer gap. The energy of the soft X-ray photons is 
determined by the backflow of the primary electrons/positrons. Each of these backflow 
particles can still maintain about 10. 6P^/^ ergs and deposit on the stellar surface. This 
energy will be emitted as soft thermal X-rays from the stellar surface (Halpcrn & Ruderman 
1993), whose characteristic energy is given by Ex{fo) ~ 1.2x lOVo^^^if p-5/i2gy_ Although 
the thermal X-ray photon density is low, every pair resulting from X-ray and high-energy 
photon interactions can emit ~ 10^ high-energy photons when they are accelerated in the 
gap. Such a large multiplicity can produce sufficient number of pairs as to sustain the 
outer gap. Assuming r — Rl/'^ and head-on collision, they obtained the fractional size of 
the outer gap as 

/„(S,P)«.5.5P^«/^^</^ (3) 

from the condition of photon-photon pair production E^E^^l — cos{9x-y)) = 2(meC^)^. 

However, this model did not take into account the fact that when the magnetic incli- 
nation angle (a) becomes large, the characteristic energy of high-energy photons from the 
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gap, which depends on a, also increases. In fact both and Oxj also depend on a. Most 
important, Zhang and Cheng (1997) has assumed a typical distance (r = ^) to represent 
the outer gap. In other words, most radiation and pair production activities are assumed 
to take place at this characteristic region. They assume that when /^(r = ^,B,P) > 1. 
the outer gap does not exist. However, even the mid-distance to the outer gap depends on 
the inclination angle because the outer gap begins at the null charge surface, which depends 
on a. Realistically, as long as pairs can be produced beyond the null surface, the outer 
gap should be still active, namely, accelerating charged particles. Of course, for regions with 
/(r) > 1. the gap potential should drop off rapidly and hence becomes unimportant. There- 
fore the active region of the outer gap should begin at null surface and stop at a distance 
Tft where /(r^) = 1. (Zhang et al. 2004). In explaining the detail 7-ray spectrum of a given 
pulsar, it is important to know the radiation coming from which parts of the pulsar mag- 
netospherc (Cheng, Rudcrman & Zhang 2000). However, in order to study the statistical 
properties of pulsars a representative region is a very useful approximation. They assume 
that the representative region of the outer gap is the averaging distance to the gap, they 
obtain the mean fractional size of the gap /(P, B, a), which can be approximately expressed 
as 

f{a,B,P)^fo(B,P)rj(a,B,P) (4) 

where r){a,B,P) is a monotonically function oi a, B and P. It roughly decreases by a factor 
3 from young pulsars with large inclination angle to old pulsars with small inclination angle. 
The integral expression of /(a, B, P) is given by Eq. (38) and the variation of r] for different 
pulsars is given in figure 3 of Zhang et al. (2004) respectively. Although the variation of 77 
is only a factor of 3, the implications are very important. First the cut-off period of 7-ray 
pulsars is about a factor of 3^^^^^ = 2.4 longer for fixed a and B in comparing with the 
old model. Secondly the cut-off age of 7-ray pulsars is about a factor of 5.8 longer as well. 
This means that there are a lot more 7-ray pulsars in high galactic latitude than previously 
expected. These 7-ray pulsars with older age can move up to the high galactic latitude and 
may contribute to the unidentified EGRET sources. 



3. Gamma-ray Emission Properties of Mature Pulseirs 

For the thick outer gap (Zhang & Cheng 1997), 7-rays are produced inside the thick 
outer gap by curvature radiation from the primary pairs along the curved magnetic field 
hues. However, Cheng & Zhang (1996) studied the radiation from the charged particles 
in the curved magnetic field, and pointed out that the radiation should be described more 
accurately by a general radiation mechanism called synchro-curvature radiation mechanism, 
in which the radiation is being emitted by the charged particles moving in a spiral trajectory 
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along the curved magnetic field lines. This mechanism differs from synchrotron and curva- 
ture mechanisms in general, but reduces to either synchrotron radiation when the radius of 
curvature of the local magnetic field lines is infinite or to curvature radiation when the pitch 
angle is zero. In fact, when the synchrotron gyro-radius = '-fmc^ sin 6p/eB{r) and the 
curvature radius of field s ~ y/r^Ri is comparable, the synchro- curvature mechanism really 
provides a significant improvement, where 7 is the Lorentz factor of the accelerated particles 
and dp is the pitch angle of the charged particles in the curved magnetic field. Zhang & 
Cheng (1997) used this mechanism to describe the production of non-thermal photons from 
the primary pairs along the curved magnetic field lines in the outer gap. The primary 
pairs have an approximate power-law distribution inside the outer gap because the energy 
and density of the primary pairs depend on local values of magnetic field, electric field 
and radius of curvature. In steady state, the energy distribution of the accelerated particles 
in the outer gap is (dN/dEe) oc Ee where E^ is the energy of the accelerated particle. 
Using the general formula of the synchro-curvature radiation power spectrum given by Cheng 
& Zhang (1996) and {dN/dEg)dE(, = {dN/dx)dx, where x = s/tRl, the differential flux at 
the Earth is (Zhang & Cheng 1997) 
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dx , (5) 



where AQ is the solid angle of 7-ray beaming, d is the distance to the pulsar, Nq — 

VSe^^oNo/hRL, No ^ 1.4 x W^'f{Bu/P)Rl 70 ~ 2 x 10V'/'(5l2/P)'/^ Rc = xR^/[{l + 
rB/{xRL))cos^9p + {RL/rB)xsm''ep], Q2 = (1 / xRlMt- b / xRl) + 1 - 3{RL/rB)x) cos^ Op + 
3{Rl / r b) X cos^ 9p + (RL/r-BYx^ sin^^Ji/2, and sin^p ^ 0.79/i/25,f /V^/^aj^^/^ F{y) = 
K^/3{z)dz, where -ft's/a is the modified Bessel functions of order 5/3, y = E^/Ef. and 
Ec = {3/2){hc'j^/x){xQ2) is the characteristic energy of the synchro-curvature photons. 
x^in and Xjnax Sive the minimum and maximum values of x. For x^in, it can estimated as 
Xmin = {fin/ RlY^"^. The maximum value of x can be estimated by Xmax ~ ^A/ry/i?^ (Arons 
& Scharlemann 1979), where 9p is the angular width of the polar cap, ^/ is the angular width 
between the magnetic axis to the upper boundary field lines of the outer gap and r/ is the 
radial distance to the upper boundary field lines intercepting with the light cylinder. Zhang 
& Cheng (1997) has discussed the possible value of x^ax- In general the 7-ray beaming 
solid angle should be different for various 7-ray pulsars, which is a function of the magnetic 
inclination angle as well as the size of the outer gap(Zhang et al. 2004). Some approximate 
forms of beaming solid angle have been given (Yadigaroglu & Romani 1995; Zhang, Zhang 
& Cheng 2000). How accurate of these approximation forms are not known. For simplicity, 
in this paper, we will treat these two parameters, AQ and Xmax as constants. By fitting the 
7-ray spectrum of known 7-ray pulsars, these two parameters are chosen to be AQ ~ Isr 
and Xjnax ~ 2 respectively (Cheng & Zhang 1998). 
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In order to compare with the EGRET data, the integral flux is necessary and is given 



where E^nax is the maximum energy of gamma-rays and is chosen to be lOOGeV. 

In Fig. 1 - Fig. 3, the model 7-ray spectra are calculated for various periods, magnetic 
fields and inclination angles respectively. We can see that the spectra become softer when the 
period decreases, the magnetic field increases or the inclination angle decreases in the energy 
range of EGRET (lOOMeV-lOGeV). Alternatively speaking the spectral break becomes larger 
when the period increases, the magnetic field decreases or the inclination angle increases. 
We can understand these trends by examining the behavior of the energy break at Eq. 
1. Typically most power radiated from the primary electrons/positrons come from regions 
with characteristic distance of order of ~ Rl- Substituting Eq. 3 and Eq. 4 into Eq. 1, 
we obtain the spectral break is oc rf/'^P^/'^^B^2^^^^ which explains the relations between 
spectral variations in Fig. 1-3. However, in subsequent sections we will find that the key 
difi^erences between 7-ray pulsars in the galactic plane and those in the high galactic latitude 
are: (l)galactic plane 7-ray pulsars are younger, shorter in period and have larger inclination 
angles and (2)they satisfy different relation between and Lgd- In fact, most high latitude 
7-ray pulsars satisfy oc Lgd which means / ~ 1 (cf. Fig. 6f and the lower panel of Fig. 
7). In this case, the spectral break is oc T^^i/iep-i/s gj^cg galactic 7-ray pulsars have larger 
inclination angles and shorter periods in comparing with the high galactic latitude 7-ray 
pulsars, so the spectrum of 7-ray pulsars in high galactic latitude has a softer spectrum than 
those in the galactic plane. In Fig. 4, we compare the typical spectra between the high and 
low galactic latitude 7-ray pulsars. 

4. Monte Carlo Simulation of 7-ray Pulsars in the Galaxy and in the Gould 



In order to consider the 7-ray luminosity and spatial evolution of pulsars in the Galaxy 
and in Gould Belt, the initial values of parameters of pulsar at birth, which include the 
initial position, velocity, period and magnetic field strength, are needed. The procedure of 
Monte Carlo method and the evolution of pulsar parameters are described in our previous 
works (Cheng & Zhang 1998; Zhang & Cheng 1999; Zhang, Zhang & Cheng 2000; Fan, 
Cheng & Manchester 2001; Zhang et al. 2004). Here, we briefiy describe basic assumptions 
for generating the 7-ray pulsars (radio-quiet and radio-loud) in the Galaxy as well as in the 
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Belt 



Gould Belt: 
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1 The pulsars are born at a rate {N^s ~ 1 per century) in the Galaxy. The age of the 

Gould Belt is estimated to be ~30 Myr old and the pulsars are born at a rate of ~ 20 
Myr-i (Grenier 2000). 

2 The Gould Belt has an ellipsoidal shaped ring with semi-major and minor axes equal 

to 500pc and 340pc respectively. The Sun is displaced from the center of Gould Belt 
about 200pc towards I — 130°(Guillout et al. 1998). On the other hand, there are 
other possible interpretations of the geometry of the Gould Belt. For examples, Olano 
(1982) and Moreno et al. (1999) have given smaller size for the Gould Belt. However, 
we will show that most 7-ray pulsars originated from the Gould Belt are mature pulsars 
with age nearly 1 million years. They already move very far away from the Gould Belt. 
The most important point is our solar system is enclosed in the Gould Belt, the exact 
dimensions of the Gould Belt is not very crucial in our problem. 

3 The initial position for each pulsar in the Galaxy is estimated from the distributions 

p^{z) = (l/2;cxp)exp(-|z|/zexp) and Pr{R) = (an/ R^^^) R exp(-R/Rcxp), where z 
is the distance from the Galactic plane, R is the distance from the Galactic center, 
Zexp = 75 pc. Or = [1 - e--^--/-^-p ( 1 + i?ma^/i?exp)] i?exp = 4.5 kpc and R^^e^ = 20 
kpc (Paczynski 1990; Sturner & Dermer 1996). But the initial position of each pulsar 
is assumed to be born uniformly inside the Gould Belt. 

4 The initial magnetic fields are distributed as a Gaussian in logB with a mean value of 12.5 

and a dispersion of 0.3. Since the majority 7-ray pulsars are younger than 3 million 
years old and the field does not decay in 10 Myr (Bhattacharya et al. 1992). So we 
ignore any field decay for these rotation-powered pulsars. 

5 The initial period is chosen to be Po=10ms and the period at time t is Given by P{t) = 

(Pq + 1.95 X 10~^^-B^t)^/^. We would like to remark that the initial period is not very 
crucial for our problem because the typical age of 7-ray pulsars is of order of million 
years old. 

6 The initial velocity of each pulsar is the vector sum of the circular rotation velocity 

at the birth location and random velocity from the supernova explosion(Paczynski 
1990). The circular velocity is determined by Galactic gravitational potential and the 
random velocity is distributed as a Maxwellian distribution with a dispersion of three 
dimensional velocity = y/S x lOOkm/s (Lorimer et al. 1997). Furthermore, the pulsar 
position at time t is determined following its motion in the Galactic gravitational 
potential. Using the equations given by Paczynski (1990) for given initial velocity, the 
orbit integrations are performed by using the 4th order Runge Kutta method with 



-9- 



variable time step (Press et al. 1992) on the variables R, Vr, z, Vz and 0. Then the 
sky position and the distance of the simulated pulsar can be calculated. 

7 The inclination angle (a) of each pulsar is chosen randomly from a uniform distribution 

(Biggs 1990). 

8 The following radio selection effects are used. The pulsar must satisfy that its radio 

flux is greater than the radio survey flux threshold and its broadened pulse width 
is less than the rotation period (e.g. Sturner & Dermer 1996). We calculate the 
400 MHz radio luminosity, L^qq, of each model pulseir using the following 
distribution given by (Ncirayan &c Ostriker 1990) Plaoo{P^ P) — 0.5A^e~^, where 
A = 3.6(log(L4oo/ < ^400 >) + 1-8), log < L400 >= 6.64 + (1/3) log(P/P=^), and L400 
is in units of mJy kpc^. The minimum detectable average flux density, 
Smini of a pulsar's radio survey is estimated by using the method of Sturner 
& Dermer (1996) (also see Cheng &; Zhang 1998). The pulsar which satisfies 
LAoo/d"^ > 'S'min IS Considered to be a radio-dctcctablc pulsar, where L400 is the radio 
luminosity at 400 MHz and d is the distance to the pulsar. The radio beaming fraction 
can be expressed as (Emmering & Chevalier 1989) /r(cf) = (1 — cosa;) + (7r/2 — a;) sin a;, 
where u; — 6°. 2 x P"^/^ (e.g. Biggs, 1990) is the half-angle of the radio emission 
cone. Then, following Emmering & Chevaher (1989), a sample pulsar with a given 
period P is chosen in one out of fr{P)~^ cases using the Monte Carlo method. There 
Eire growing evidence that very strong surface magnetic field exists on the 
surface of neutron stars. For examples, absorption/emission line features 
have been observed in isolated pulsars/neutron stars, e.g. IE 1207.4-5209 
(Sanwal et al. 2002), PSR 1821-24 (Becker et al. 2003) and RBS1223 
(Haberl et al. 2003). These line features imply that the surface magnetic 
fields cire one to two order of magnitude higher than that of the dipolcir 
field. In fact, Cheng &c Zhang (1999) have analyzed the X-ray emission 
from the polcir cap regions of the rotation-powered pulseirs and concluded 
that there exists a chciracteristic surface field with strength ~ lO^^G. Such 
strong surface magnetic field can easily ciffect the beaming direction of radio 
wave. Therefore, we assume that the beaming direction of 7-rays is independent of 
the beaming direction of the radio. In general, it is possible that there is only one 
beaming, either the radio or the gamma-ray beam, pointing toward us. In fact this is 
one of key differences between outer gap model and polar gap model. The solid angle 
of 7-rays is taken to be Isr. 

9 The 7-ray threshold varies over the sky. Yadigaroglu & Romani (1995) used a fiux 

threshold of 3 x 10~^° erg cm~^s~^, which can compare to the faintest 5a sources in 
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the first EGRET catalog(Fichtel et al. 1994). However, in the third EGRET catalog, 
the faintest source in the catalog with significance \/TS >4 has a photon fiux of 
(6.2 ± 1.7) X 10"^ cm-^s-i (Hartman et al. 1999). Gonthier et al. (2002) have argued 
that this threshold could be reduced for | & |> 10°. In our analysis, we include the 
criterion of the likelihood \/TS >5 (~ 5a) which corresponds to the energy threshold 
of S^^{E^ > lOOMeV) > 1.2 x 10"^° erg cm-^s"! for | 6 |< 10° but decreases the 
threshold to Stj'iE^ > lOOMeV) > 7.0 x 10"^^ erg cm-^g-i for | 6 |> 10° (Gonthier et 
al. 2002). 



5. Simulation results and Discussion 

Following the procedure described in section 4, we perform Monte Carlo 
simulation of Galactic pulsars born during past 3 x 10"^ yrs. We use the code of 
Cheng &c Zhang (1998) (also see Zhang et al. 2000) in our simulation, in which 
Pcirkes Multibeam survey is not included. We summarize the various components of 
simulated 7-ray pulsars in table 1. Our results suggest that there should be more unidentified 
EGRET 7-ray sources identified as radio pulsars. In particular, the Parkes multi-beam pulsar 
survey is expected to discovery more radio pulsars (Manchester et al. 2001). Torres et al. 
(2001) have found five new possible radio pulsar-Unidentified source associations. Most 
recently, Kramer et al. (2003) have used the newly release survey results to correlate the 
unidentified EGRET 7-ray sources and, have found more than 35 coincidences and around 
20 probable associations. It is important to note that the number in each box of table 1 is 
sensitive to the input parameters, e.g. detection threshold, birth rate, initial distributions 
of pulsars etc., for a given theoretical model. However, their ratios are less sensitive to most 
input parameters except their relative birth rate. It is interested to note that the 7-ray 
pulsars in the Gould Belt contribute significant fraction of total 7-ray pulsars for | 6 |> 5° 
and become unimportant for | 6 |< 5°. This is simply because the solar system is enclosed 
by the Gould Belt. 

In Fig. 5a-e, we compare our simulated distributions of period, period derivative, mag- 
netic field, distance and energy flux with the observed data of 8 known radio-loud 7-ray 
pulsars (i.e. Crab, Vela, Geminga, PSR 0656+14, PSR 1046-58, PSR 1055-52, PSR 1509-58, 
PSR 1706-44, PSR 1951+32 ) using KS test. Although PSR 0656+14 and PSR 1509-58 are 
confirmed in 7-ray band, there are insufficient 7-ray photons to provide the information of 
pulsed fraction. So we did not include these pulsars in the cumulative plot of energy fiux. 
The maximum deviations of period, period derivative, magnetic field, distance and 7-ray 
energy flux distributions from the observed distributions are 0.36, 0.32, 0.25, 0.33 and 0.21 
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respectively. It can be seen that four of five accumulative distributions cannot be rejected at 
better 80% confidence level, and period cumulative distribution cannot be rejected at better 
than 90% confidence level. Therefore we conclude that the model results do not conflict with 
the observed data of 7-ray pulsars. 

In Fig. 6a-f, we plot the normalized distributions of period, magnetic field, age, distance, 
inclination angle and the fractional size of outer gap. The solid lines and the dashed lines 
are 7-ray pulsars located at | 6 |< 5° and | 6 |> 5° respectively. Wc can clearly see that 
there arc two classes of 7-ray pulsars. 7-ray pulsars located at | 6 |< 5° have shorter periods, 
younger, larger inclination angles and smaller outer gap size in comparing with 7-ray pulsars 
located at | 6 |> 5°. As we have mentioned that 7-ray pulsars in high galactic latitude are 
dominated by mature pulsars, which are old enough to move up the high latitude and to 
evolve to longer periods. Young 7-ray pulsars will be dominated in the galactic plane, which 
have stronger 7-ray luminosity, larger inclination angle, shorter period and larger magnetic 
field. In Fig. 6, 7-ray pulsars in high galactic latitude are actually closer than those in 
galactic plane. It is because they are older and weaker 7-ray pulsars, so they must be near 
otherwise they cannot be detected. 

In Fig. 7, we plot versus Lsd, where the upper panel is 7-ray pulsars located at 

I 6 |< 5° and lower panel is 7-ray pulsars located at | h |> 5°. Wc can see that there are 
two distinctive regions. For Lgd > 3 x 10'^'^ erg/s, the relation between L^, and Lsd is rather 
scattered. However, for Lgd < 3 x 10^^ crg/s, is proportional to Lgd- In particular, the 
7-ray pulsars in high latitude, most pulsars satisfy oc L^^, which means / ~ 1. This is 
supported by Fig. 6f, in which the fractional sizes of outer gap / for 7-ray pulsars located 
at I 6 |> 5° are all close to unity whereas the fractional sizes of outer gap for 7-ray pulsars 
located at | 6 | < 5° have a wide distribution. If oc L^^ is used to fit all 7-ray pulsars, 
then 13 ~ 0.3 for 7-ray pulsars located at | 6 |< 5° and (3 ~ 0.6 for 7-ray pulsars located at 
I & |> 5° respectively. Since the outer gap size is a function of -B, Panda. For young pulsars, 
there are more combination of these three parameters to make the outer gap size less than 
unity. On the other hand, for old 7-ray pulsars their periods are already so long that there 
are not much room for other two parameters to vary the outer gap size away from unity. 
The spatial distributions of 7-ray pulsars arc given in Fig. 8. 

In Fig. 9a-d, we compare distributions of 7-ray pulsars from the Galaxy (dotted lines) 
and from the Gould Belt (solid lines). We can see that the properties of 7-ray pulsars from 
two sources are very similar except 7-ray pulsars from the Gould Belt have lower population 
and closer. These two properties results from the lower birth rate in the Gould Belt and the 
solar system is inside the Gould Belt. 



6. Conclusion and discussion 



We have studied the 7-ray emission properties of pulsars in the Galaxy as well as in 
the Gould Belt by using a new self-consistent outer gap model, which includes the effects of 
inclination angle and average properties of outer gap. We have found that this new model 
can produce more mature 7-ray pulsars, which have longer period, smaller inclination angle 
and typical age of one million years old, than the old model. In fact the mature 7-ray pulsars 
dominate in the high galactic latitude. The spectra of these mature pulsars are significantly 
softer and weaker than the young pulsars in the galactic plane. It is because the spectral 
hardness is roughly determined by the typical photon energy given in equation (1), which 
characterizes the position of spectral break. The explicit form of how the typical photon 
energy depends on the inclination angle is given in equations (28) and (38) of Zhang et 
al. (2004). Although the form is quite complicated, we can roughly understand why smaller 
inclination angle gives smaller typical photon energy and hence softer spectrum by equation 
(1). It is because the typical position of the outer gap is the null surface where becomes 
larger for smaller inclination angle. Prom the simulation results, it is not clear why 7-ray 
pulsars in high latitude should have smaller inchnation angle. It is because all of these high 
latitude pulsars are older pulsars, which have longer period. Without the compensation of 
the effect of inclination angle, the fractional size of the outer gap / is larger than unity ( 
cf. equation 3). In order to maintain the outer gap size less than unity, rj in equation (4), 
which contains the effect of inclination angle must be small. In fact rj becomes small for 
small inclination angle (cf. Fig. 1 of Zhang et ai. 2004). Therefore high latitude gamma-ray 
pulsars tend to have smaller inclination angle and harder spectrum. On the other, this seems 
to imply that the confirmed radio-loud 7-ray pulsars in the galactic plane should have larger 
inclination. This statistical correlation has not yet been established. But we would hke to 
point out that the values of inclination angle of pulsars are very difficult to be measured 
accurately and there are only 7 confirmed 7-ray pulsars. So it is very difficult to find such 
correlation in radio-loud 7-ray pulsars. However, we do believe that this correlation should 
be there when more radio-loud 7-ray pulsars are discovered in the galactic plane. 

We have also used a Monte Carlo method to simulate the properties of the 7-ray pulsar 
population in the Galaxy as well as in the nearby Gould Belt in terms of the revised outer 
gap models. The initial magnetic field, spatial and velocity distributions of the neutron star 
at birth which are obtained by the radio pulsar statistical studies have been used in our 
simulations. We have obtained the spatial, distance, period, age, magnetic field, inclination 
angle, photon flux distributions of the radio-loud and radio-quiet 7-ray pulsars. We flnd 
that the properties of 7-ray pulsars between | & |< 5° and | 6 |> 5° are very much different. 
Galactic plane 7-ray pulsars are younger, shorter in period and have larger inclination angles. 
They satisfy different relation between and Lsd- The high latitude 7-ray pulsars satisfy 
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oc Vlf but the galactic plane pulsars satisfy oc U^f. 

The present model predicts very similar numbers of 7-ray pulsars as old models but many 
more high latitude 7-ray pulsars than old models (Cheng & Zhang 1998; Zhang, Zhang & 
Cheng 2000). It is because the new model allows the outer gap can exist about a few miUion 
years old for appropriate combination of B,P and a, so old 7-ray pulsars can move up to high 
galactic latitude. Furthermore, the nearby Gould Belt also contribute significant number of 
7-ray pulsars in high latitude. Torres et al. (2003) have predicted that AGILE can detect 
more 7-ray pulsars (for a general review of AGILE cf. Tavani et al. 2001). Perhaps AGILE 
can provide some clues to differentiate polar gap model predictions (Gonthier et al. 2002) 
and outer gap model predictions, and GLAST makes the final verdict. Finally we would 
like to remark that the better predictions on how many 7-ray pulsars should be detected by 
AGILE and GLAST should include a better expression of 7-ray sohd angle, which depends 
on the inclination angle as well as the outer gap size. However, in order to have a reliable 
expression of 7-ray solid angle, we need to carry out three-dimensional model calculation 
and come up with an approximate expression from Monte Carlo simulation. 

We thank an anonymous referee for the very useful comments. This work is partially 
supported by a RGC grant of Hong Kong Government, 'Hundred Talents Program of CAS' 
and the National 973 Projection of China (NKBRSFG 19990754). 
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Table 1: Various Components of Simulated 7-ray Pulsars 





Gould Belt 


Galaxy 


Birth rate 


l/50,000yr. 


1/lOOyr. 


Galactic latitude b 


< 5° 


> 5° 


< 5° 


> 5° 


Radio-loud 


1 
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12 


4 


Radio-quiet 


2 


13 


40 


22 



10' 10- 

E (MeV) 



10= 



Fig. 1 Distribution of the energy flux vs energy of the photon. In this flgure, Bu has been 
taken to be 3, P to be 0.4 s and a to be 20° (sohd hne), 40° (dashed hne) and 60° (dotted 
hne) . 

10-- E ' ' ' j 



10= r 




E (MeV) 

Fig. 2 Distribution of the energy flux vs energy of the photon. In this figure, P has taken to 
be 0.4 s, a to be 40° and B12 to be 1 (sohd line), 3 (dashed line) and 5 (dotted line). 
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Fig. 3 Distribution of the energy flux vs energy of the photon. In this figure, B12 has taken 
to be 3, a to be 40° and P to be 0.2 s (sohd fine), 0.4 s (dashed fine) and 0.6 s (dotted fine). 



10" 




Fig. 4 Distribution of the energy flux vs energy of the photon. The sohd hne is the typical 
photon spectrum of a galactic 7-ray pulsar with B12 = 5, P = 0.1s and a = 70°, and the 
dashed line is the typical photon spectrum of a 7-ray pulsar in high galactic latitude with 
B12 — 2, P — 0.3s and a — 20° respectively. 
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Fig. 5 Normalized cumulative distributions (dashed curves) of period, period derivative, 
magnetic field, distance and 7-ray energy fiux of 7-ray pulsar population for our model with 
radio selection effects. For comparison, corresponding distributions (sohd hnes) of the 
observed 7-ray pulsars are also shown. 
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Fig. 6 Normalized cumulative distributions of period, magnetic field, distance, age, inclina- 
tion angle and the averaged fractional size of outer gap of both simulated 7-ray pulsars in 
|6| < 5° (solid lines )and in \b\ > 5° (dashed lines), which are labelled by A, B, C, D, 
E, F respectively. 
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Fig. 7 Plot of 7-ray luminosity vs. spin-down luminosity for the simulated 7-ray pulsars. 
Upper panel for the case in |6| < 5°, and bottom panel for the case in |6| > 5°. 
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Fig. 8 Normalized distributions of the simulated 7-ray pulsars in Galactic longitude and 
latitude. 
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Fig. 9 Normalized distributions of period, magnetic field, distance and age of the simulated 
7-ray pulsars for |6| > 5°. Dotted histograms represent the distributions of the simulated 
7-ray pulsars without those produced in the Gould belt, and solid histograms represent the 
simulated 7-ray pulsars only produced in the Gould belt. 



